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Abstract 

The structural, dielectric and dynamical properties of the low temperature antiferromagnetic 
orthorhombic phase of CaMnOa have been computed from first principles using a density functional 
theory approach within the local spin density approximation. The theoretical structural parameters 
are in good agreement with experiment. The full set of infrared and Raman zone-center phonons 
is reported and compared to experimental data. It is shown that coherently with the anomalous 
Born effective charges and the presence of low frequency polar modes, the static dielectric constant 
is very large and highly anisotropic. 
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I. INTRODUCTION 



Magnetic perovskite oxides constitute an exciting subject of study. Indeed, the interplay 
between their structural, magnetic and transport properties make them fascinating from 
both experimental and theoretical point of views. For example, the mixed valent perovskite 
Cai-xLuxMnOs is one of the most studied material for its colossal magnetoresistance[l|. 
Also, the oxygen deficient manganites LaMnO^s and CaMnO^^s are interesting for their 
transport and optical propertiesp]. 

In this paper we report a first-principles study of stoichiometric CaMnO^. The ground- 

r 

state crystal structure of CaMnO^ is orthorhombic with space group Pnma[4]. The structure 



bur formula per unit 



im- It 



can be regarded as a distorted perovskite structure having 

is an insulator with an observed bandgap of about 3 evjoj. The magnetic structure is 
antiferromagnetic and the G— type order is energetically the most favourable. The magnetic 
interactions between the Mn ions are due to superexchange interactions. Because of Mn~^^ 
configuration, there is no tendency to Jahn- Teller distortion. The observed Neel temperature 
is about 130 K, which estimates the exchange energy to be 6.6 meV. 

Previous first-principles calculations on CaMnO^ focused mainly on the electronic struc- 
ture of the system. In this paper we extend the study and report results concerning the 
dielectric and dynamical properties of is ground-state orthorhombic phase. 

II. TECHNICAL DETAILS 

The first-principles simulations were performed according to the density functional theory 
scheme (DFT) within the local spin density approximation (LDA) and using the plane-wave 



implementation of the ABINIT package [9j]. We used the Hartwigsen-Goedecker-Hutter [11 1 
(HGH) parametrization for the pseudopotentials, where the 3s and 3p orbitals were treated 
as valence for Mn and Ca atoms and 2s and 2p orbitals were considered as valence for O 
atoms. The total number of valence electrons is therefore 15 for Mn, 10 for Ca and 6 for 
oxygen. Convergency was reached for an energy cutoff of 72 Hartree for the plane-wave 
expansion and a 6 x 4 x 6 A;-point mesh for the Brillouin zone integration. The phonon 
frequencies. Born effective charges and dielectric tensor were computed according to the 
density functional perturbation theory ^10] (DFPT) scheme as implemented in the ABINIT 
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package. 



III. STRUCTURAL PROPERTIES 

At high temperature, CaMnOs adopts a cubic perovskite structure. At room tempera- 
ture, it crystahzes in an orthorhombic Pnma (N.62) structure with 20 atoms in the primitive 
unit cell and adopts a G— type antiferromagnetic configuration [4]. This orthorhombic phase 
results from the relaxation of three kinds of antiferrodistortive (AFD) instabilities which 
are related to rotations of oxygen octahedras around the Mn-0 axis. In Glazer's notation, 
these tilts are described by a^b+c^. Since these oxygen rotations do not strongly affect 
the high symmetric structure, the orthorhombic phase of CaMnOs can be described as a 
pseudocubic structure, the volume of which can be estimated as v^acX2ac x v^Sc where a.^ 
is the pseudocubic cell parameter. 

In our calculations, we imposed a G— type antiferromagnetic order and did the structural 
relaxation at fixed volume. The volume used corresponds a pseudocubic cell parameter of 
3.73A which corresponds to the experimental cubic lattice constant and also closely agree 
with the experimental pseudo-cubic lattice constant of the orthorhombic phase ^. The cell 
shape and atomic positions were relaxed according to this constraint and the results are 
reported in TABlB 

The relaxed structure is in excellent agreement with the experimental data. The cell 
parameter a and h are only slightly overestimated while the c parameter is slightly under- 
estimated (errors of ~ 0.5 %). The calculated atomic positions are also in good agreement 
with the experimental data, the distortion with respect to the ideal cubic positions being pre- 
dicted with an accuracy similar to what is usually achieved in the class of ABO3 compounds 

H. 

From the inspection of the density of states, we found a theoretical electronic bandgap 
of 0.8 eV for our relaxed antiferromagnetic orthorhombic structure. Although this strongly 
underestimates the experimental value (3.1 eV [ifil), as usual within the LDA, our calculation 
properly reproduces the insulating character of the structure. 
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TABLE I: Theoretical and experimental unit-cell parameters (A) and non-equivalent atomic po- 
sitions (reduced coordinates) in the antiferromagnetic orthorhombic Pnma structure of CaMnOs. 
The Wyckoff position of the reported atoms are mentioned in brackets. 

IV. DIELECTRIC PROPERTIES 

The Born effective charge tensors have been calculated for the the four atoms of TAB. [H 
The full tensors are as follows : 
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Z*{Ca) 



( 2.51 0.00 0.19^ 
0.00 2.43 0.00 
^0.27 0.00 2.52^ 
[2.29 2.74 2.43] 



Z*[Mn) 



( 6.82 -0.37 -0.72 ^ 
0.07 5.78 1.42 
-0.69 -1.54 6.56 
[7.41 6.05 5.70] 



^^-1.73 0.00 0.06 ^ 

0.00 -5.11 0.00 
^-0.20 0.00 -1.76 / 
1.67 - 1.82 -5.11] 



( -Zm -0.08 -1.94^ 
-0.11 -1.55 -0.01 
^-1.97 0.00 -3.61 
5.66 - 1.77 - 1.53] 



I 



Below each tensor, the main values of the symmetric part of Z* is also mentioned, for 
direct comparison with (i) the nominal atomic charges (ZAr(Ca) = +2, Zf^iMn) = +4 and 
Zn{0) = —2) and (ii) the Born effective charges in the cubic phase at = 3.73A(Z^(Ca) = 



'he trends are rather 
1^ . As usual within 



+2.61, Z^{Mn) = +7.43, Z^{0^) = -2.55 and Z^{0\\) = -4.94). 
similar to what was previously re por ted for the isostructral CaTiO^ 

this class of ABO3 compounds 18|, ll9|] , we notice the anomalously large values of the Mn and 
O charges. As usual also, in comparison to the highly symmetric cubic structure, there is a 
reduction of the anomalous values in the low symmetry phase, at the exception of Z*y{Oi). 
Finally, we notice the existence of small asymmetric contributions, as allowed by symmetry 
in the orthorhombic phase. 



Although its value is not expected to be very accurately predicted within the LDA [20[, 
we also report the calculated optical dielectric tensor : 

/ 11.3 0.0 0.0 ^ 
0.0 13.1 0.0 
0.0 0.0 10.8 



V 
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The calculated tensor is only slightly anisotropic. Its elements are significantly larger 



than in related compounds like CaTiO^ (eoo ~ 6 [IJ]), coherently with the smaller LDA 
bandgap of CaMnOs. The optical dielectric tensor in the orthorhombic phase does not 
differ significantly from that previously reported in the ideal cubic structure (eoo = 10.43 



from our calculations and 



11.25 from Ref. 



21I ). We notice however that it is slightly 



larger than in the cubic phase, which is a trend different from that reported for orthorhombic 



CaTiO^ in which decreases when non-polar distortions are frozen into the structure 



1^ 



V. DYNAMICAL PROPERTIES 



The irreducible representation in the orthorhombic Pnma phase CaMnOa at the F point 

is: 

7Ag © ^Big © 752,; © 553,; © 10Si„ © S^a^ © lO^s^ © 8A„ 

Over this decomposition, 3 modes are acoustic (symmetries i?2u, -Bsm), 8 are silent 
(symmetry y4„), 24 are Raman (R) active (symmetries Ag, Big, and B^g) and the 
last 25 modes are infrared (IR) active (symmetries i?2« and i?3„). According to the 
structure defined in TABJH the IR B^u modes are polarized along the x direction, B2u along 
y direction and Bi^ along z direction. The eight silent modes, not further discussed 
below, are calculated at the frequencies of 123, 140, 179, 220, 313, 392, 438 and 466 cm~^ 

A. Raman active modes 

We report in TAB Hi] the calculated frequencies of the Raman active modes and a re- 
lated new assignment of the experimental data. The latter is compared to the assignment 



previously proposed from shell-model results in Ref. 15|. 

Experimentally, the main observed modes are those of Ag symmetry. 10 A^ modes were 
measured at 150, 160, 184, 243, 278, 322, 382, 438, 487 and 615 cm'^ while only 7 modes 
are expected from group theory. The frequency at 615 cm^^ was unambigously assigned to 
impurity and the modes at 382 and 438 cm^^ were kept aside since their frequencies differ 
strongly from those of the shell model. Below 200 cm^^ the assignment was particularly 
ambiguous since two modes were calculated at 154 and 200 cm^^, while three were oberved 
in the spectra. 

From our calculation, the experimental lines at 243, 278 and 322 cm~^ can be assigned to 
the Ag modes at 250, 275 and 314 cm~^ respectively. For the low frequencies, the assignment 
can be significantly improve if we assign to the calculated modes at 152 and 167 cm~^, 
the experimental modes at 150 and 160 cm~^ rather than those at 160 and 184 cm~^ as 
previously proposed (i.e. rulling out the experimental mode at 184 cm"^ instead of the one 
at 150 cm^^). For the high frequency A^ modes, the overall agreement between experimental 
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and theoretical data is the best if we assign the experimental line at 438 cm~^ (not assigned 
in Ref ■ [isl ) to the calculated mode at 450 cm~^ and the one at 487 cm~^ to that calculated 
mode at 504 cm~^. 

The experimental line at 179 cm^^, assigned to a Big mode in Ref. [isl . is in good agree- 
ment with our calculation (189 cm^^). Amongst the two experimental lines at 320 and 
564 cm~^ previously assigned to of B^g modes 15||, the first one is reproduced within our 
calculation with a very good accuracy (320 cm~^), contrary to the second for which we get 
a theoretical frequency of 469 cm~^). Since the assignment between Big and B^g modes is 
rather ambiguous, we propose that the experimental line at 564 cm~^ might in fact corre- 
spond to the Big mode calculated at 595 cm^^. 



Finally, from Ref. 



15|, it seems that the B2C, modes are theoretically predicted with less 



accuracy since modes calculated with the shell model at 180 and 425 cm~^ are associated to 
the lines measured at 258 and 465 cm^'^. From our first-principles calculation, this agreement 
seems significantly better and we propose to assign the experimental line at 258 cm~^ to the 
theoretical mode at 227 cm~^ and the line at 465 cm~^ to the theoretical mode at 465 cm~^. 



B. IR active modes 

Table IIIII summarizes the calculated IR frequencies of the TO modes. Unfortunately 
no IR measurement on monocrystal in the orthorhombic phase of CaMnOs was found to 



compare with our results. An IR spectra was reported in Ref.[l7j but on a polycrystalline 
sample and no symmetry attribution was reported for the 15 observed frequencies, which 
makes comparison very difficult considering the high number of modes. The only remark 
that can be done concern the high frequency part of the spectrum. The three highest 
frequencies are measured in Ref.[17|] at 533, 580 and 628 cm~^, which deviate strongly with 
our results where the maximum frequency of IR modes is calculated at 504 cm~^. This could 
eventually be related to inaccuracies in our calculations but our calculation also suggest that 
these experimental frequencies might also likely correspond to combination of modes (for 
exemple, the highest measured frequency 628 cm~^ can be recover as being the exact sum 
of the modes calculated at 305 cm~^ and 323 cm^^). 



In addition to the frequencies, we also report in TAB IIIII the mode effective charges ( 22|). 



the oscillator strengths and the contribution of each polar mode to the static dielectric tensor. 
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Symmetry Present Reference [15j 
FP Exp. [15] SM Exp. 
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TABLE II: Comparison of the calculated and experimental frequencies (cm^^) of the Raman modes 
of the antiferromagnetic orthorhombic phase of CaMnO^. The first column labels the symmetry 
of each mode. The second and third columns correspond to the present first-principles (FP) 
calculations and the related assignment of the experimental data of Ref . UM ■ The fourth and fifth 

n 

columns corresponds to the shell-model (SM) calculations and related assignment in Ref.|15l]. 
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Symmetry 
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TABLE III: Symmetry and frequency of the 25 IR modes. The mode effective charges (Z^^,^, as 



defined in Ref. [2^) and the oscillator strengths (/S^") are also provided for each mode as well 

as the contribution to the static dielectric constant egm- From our conventions, aa = xx for B^u 

modes, aa = yy for i?2u niodes and aa = zz for Biu modes. 
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The total static dielectric tensor can be decomposed as follows: 




(1) 



m 



where a and (3 are the cartesian directions (x, y or z), eoo is the optical dielectric tensor and 
eo,m is the contribution to the dielectric constant of each individual phonon mode m. This 
latter contribution is computed from the following relation: 



where Vt is the volume of the cell, S'^ and Um are respectively the oscillator strength and 
the frequency of the mode m. 

The computed static dielectric tensor is reported below : 



Not only the static dielectric tensor takes values significantly larger than the optical 
dielectric tensor but, contrary to the latter, it is also highly anisotropic. Some insight into 
this result is provided from the inspection of TAB lllli Along the x direction, the static 
dielectric constant [e^^ = 262) appears as mainly due to the contribution of the lowest B^u 
mode at lOlcm"^ (238), which combine giant mode effective charge and oscillator strength 
and a low frequency. Along the z direction, there are still highly polar mode (Bi^) with 
giant mode effective charges but at higher frequencies so that the static dielectric constant 
is smaller. Along the y direction, the most polar mode (B2u) is still at larger frequency. 

VI. CONCLUSION 

The structural, dielectric and dynamical properties of the anti-ferromagnetic orthorhom- 
bic phase of CaMnO^ have been studied from first-principles. The relaxed structure is 
in good agreement with experimental data. As within the cubic phase, the Born effective 
charges of Mn and O are highly anomalous. The whole set of zone-center phonon modes has 
been computed and a new assignment of experimental data as been proposed. The static 
dielectric tensor has also been obtained : it shows amplitudes comparable to CaTiOs and 




(2) 



/ 262 0.0 0.0 \ 



Co = 



0.0 68 0.0 



y 0.0 0.0 120 ) 
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is highly anisotropic. Both the inspection of the phonon modes and static dielectric tensor 
emphasize that, although from the structural point of view the orthorhombic phase can be 
considered as a pseudo-cubic structure, from the dielectric and dynamical point of view, it 
is highly anisotropic. 
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